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Abstract. The momentum of an MHD wave has been examined from the view 
point of the electromagnetic momentum expression derived by Minkowski. Basic 
calculations show that the Minkowski momentum is the sum of electromagnetic 
momentum and the momentum of the medium, as proposed in some of the past 
literature. The result has been explicitly confirmed by an example of an MHD wave, 
whose dynamics can be easily and precisely calculated from basic equations. The 
example of MHD wave also demonstrates the possiblility to construct a symmetric 
energy-momentum tensor based on the Minkowski momentum. 
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1. Introduction 

The Minkowski- Abraham controversy has been discussed by a number of authors 
over a hundred years. Minkowski [IJ proposed the electromagnetic momentum 
density in a dielectric medium must be D x B, and Abraham 121 [3l proposed E x H 
for that (in the present letter symbols have conventional meanings, e.g., E = electric 
field, unless otherwise stated). There have been published numerous papers on this 
problem both theoretically and experimentally, but the final conclusion is still yet to 
come; papers are still beeing published in this century (see, e.g., ||4| for a review). 

Several authors BHUH, HI pointed that the electromagnetic field inevitably affect 
the dynamics of the medium to change its energy-momentum, and therefore, the 
energy-momentum of an electromagnetic wave must include the contribution of the 
medium. In the present letter we show that the Minkowski momentum is the sum of 
electromagnetic momentum and the momentum of the medium. Feigel lUl obtained 
a similar result based on the Noether's theorem using Lagrangian formulation. 
Compared to his elegant approach, the calculation here is rather a down-to-eath 
type, which is more closer to the Minkowski's original derivation. This approach is 
less elegant, however, easier to understand its meaning intuitively. 

Perhaps the largest weak point of the Minkowski momentum is the fact that the 
four dimensional energy-momentum tensor does not become symmetric with this 
momentum, which means the violation of angular momentum conservation (see, e.g. 
0). Most of the past literature argued the legitimacy of the momentum part of the 
tensor in this point. Here, in contrast, we elucidate the possibility to alter the energy 
part to make a symmetric tensor; provided the momentum part of the Minkowski 
energy-momentum tensor includes the momentum of medium, the same should be 
true for the energy part. To treat it in a relativistically consistent way, the mass flux 
must be included in the energy flux even in the non-relativistic regime. The energy- 
momentum tensor with Minkowski momentum can become symmetric when the 
mass flux is taken into account. 

The consideration stated above is confirmed by an example of an MHD wave in 
a collisionless magnetized plasma. Usually the behavior of an ordinary medium 
is complicated and need to calculate microscopic states of molecules, which is 
difficult to solve exactly. A collisionless plasma is, in contrast, easy to calculate 
its response to the electromagnetic field from the classical basic equations (Maxwell 
equations and Newtonian mechanics). Here in this short letter we use the MHD 
approximation, however, if one wishes it is possible to derive an exact solution of the 
basic equation system to confirm the result. The result agree with the "frozen-in" of 
a magnetized plasma, which has been confirmed by a wide variety of experimental 
and observational facts. 
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2. Basics 



Microscopic Ampere's equation in a medium is 

-eo^+f'o'VxB=J (1) 

Suppose there is no external current, and J consists of the polarization current Jp 
and magnetization current Jm, which are generated in response to the electric field 
E and magnetic field B respectively. We introduce the polarization vector P and 
magnetization vector M such that 

^P=Jp, VxM = Jm (2) 

In this context, P and M should be understood as convenient mathematical 
expressions to represent the response of the medium to the electromagnetic field, 
rather than real physical entities. 

When averaged over a microscopically large but macroscopically small volume, 
P and M are assumed to have simple linear relations to the electromagnetic fields as 

P = XpE, M = XmB. (3) 

The linear coefficients Xp and Xm are matrices in general because the medium may 
not be isotropic (as in our example of magnetized plasmas). The fields D and H are 
then defined as macroscopic quantities as 

D = eoE + P = eE, H = fi-^B + M = fr^B (4) 



3. Minkowski Momentum 



The momentum of microscopic electromagnetic field is CqE X B and its conservation 
law is 

^(eoE X B) + V • r + Gp + Jm) X B + (Qp + Qm)E = , (5) 

where T is the Maxwell stress tensor and we denote dTij/dxi = (V • T)j in 
short. The polarization/magnetization charge Qp and Qm are the result of 
polarization/magnetization current {dQp^M/dt = VJpm)- The charge due to 
magnetization current vanishes when averaged, Qm = since ]m satisfies 

The third and fourth term of (O are the Lorentz and Coulomb force acting on the 
medium, and therefore, it can expressed by the momentum change of the medium. 

(Jp + Jm) X B + (Qp + Qm)E = + V • Tm , (6) 

where g and Tm are the momentum density and stress tensor of the medium. It should 
be noted that the right hand side of the above expression has mathematical ambiguity. 
If we define new values of momentum/stress by g' = g+aandT' = T+G with arbitrary 
vector a and tensor G that satisfy da/dt = VG = 0, they also satisfy the above equation. 
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Therefore, g and T do not necessarily have to be the total momentum/stress of the 
medium. For example, the medium may contain a part that does not interact with 
the electromagnetic field, and such part causes this ambiguity 
From dU we obtain 



Jp X B = -(P X B) + xp 

1. 



(E • V)E + -VE^ + E(VE) 



(7) 



(BV)B + -VB' 



Jm X B - Xm 

Here we neglected cross terms of fluctuation in averaging as usually done in this kind 
of calculation, e.g., P x B = P x B. Combining (H), ^ and we obtain 

0, 



^(£oE X B + g) + V • (T + Tm) = ^(D X B) + V • r 



(8) 



(9) 



where T' is the stress tensor in a dielectric medium defined as 

f ;. = EiDj + ^-%Hj - ^6ij{E-D + B • H) , 

which is the sum of the fluxes of electromagnetic momentum and momentum carried 
by the medium. 

Now that we understand the Minkowski momentum includes the part of the 
medium then so should be for the energy. The energy is equivalent to mass in 
relativity, thus the energy flux of the medium must include mass, and then the flux 
may take the form of D x B to make the energy momentum symmetric. We will check 
this for the case of an Alfven wave in the following. 



4, MHD wave 



Let us confirm the above discussion with an example of an MHD wave. Suppose a 
linearlypolarized one dimensional ((9/(9x = d/dy = 0) MHD wave (Alfven wave in this 
case) propagating in the z direction, which is the direction of the background magnetic 
field: Bq = Bz- The wave amplitude is small enough for linear approximation, and the 
plasma velocity is low enough for non-relativistic approximation. Also we assume 
"cold plasma limit", which means the thermal energy of plasma particles is negligibly 
small. 

The wave has an electric field perpendicular to its propagation direction, and we 
take the x axis in this electric field direction. Then the current is also in the x direction, 
whereas the magnetic perturbation and the plasma velocity is in the y direction (see 
Appendix). The magnetic field created by cyclotron motion of plasma particles is 
negligible in a cold plasma limit, and thus we treat the case with fj. - [Iq hereafter. 

The current is the polarization current due to the temporal change of the electric 
field, which is 

Ua d 
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where is th Alfven speed defined by Va = Bq/ ^/JioP with p being the mass density 
of the plasma. From (O and the above expression we obtain 

D,=Eo(i+^y, (11) 

The plasma frozen-in condition E + v X B = means that the plasma is moving in the 
y direction with the E x B drift speed as 

Then the momentum carried by the plasma particles is 



- 



pvy = —E, . (13) 



A 



The y component of the Minkowski momentum can be calculated from (|TT) and (|13|) , 
which is 

(D X B)y = D,Bo = eoExBo + pVy . (14) 

The above expression means the Minkowski momentum is the sum of electromagnetic 
momentum and momentum of the plasma particles as long as the frozen-in condition 
is satisfied. 

When multiplied by c^, the first term of the right hand side of ((T4)l becomes 



the electromagnetic energy flux in the y direction. The second term becomes the 
relativistic energy comes from the rest mass, which is the predominant energy flux in 
the non-relativistic limit here; thermal or kinetic energy flux is negligible. The balance 
equation of the energy-momentum tensor is in a derivative form, and therefore, it 
has an ambiguity as discussed below for the momentum. The same is true for the 
energy, and its conservation also holds when we add the mass density and mass flux 
(p, pv), since dp/dt + V(pv) = 0. The energy momentum tensor becomes symmetric 
with these terms. 



5. Concluding Remarks 

Momentum carried by an electromagnetic wave has been examined with an example 
of an MHD wave. It has been shown that the total momentum (electromagnetic 
momentum plus momentum of the medium) is expressed by D X B as proposed by 
Minkowski. This result is based on the very simple and basic two properties of an 
MHD plasma, the frozen-in condition ((T2)l and polarization current (|TOl) , namely. It 
would not be exaggeration if one says the whole kingdom of MHD plasma physics 
would fall if these two basic properties were wrong. 

Here in this letter we examined a simplest case of an parallel (to the B field) 
propagating MHD wave, but similar calculations can be done for more complicated 
plasma waves to confirm the result here. A coUisionless plasma contains a wide 
variety of wave phenomena, and the properties of waves can be precisely calculated 
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at least in the linear limit. Calculation of the Minkowski momentum for various 
plasma waves would be a good exercise to understand the Abraham-Minkowski 
controversy. 

The drawback of the Minkowski momentum has been believed that the 
momentum fails to form a symmetric four dimensional energy-momentum tensor 
when coupled with the Poyinting flux; an asymmetric energy-momentum tensor 
means the violation of angular momentum conservation. This difficulty can be 
overcome when we include the mass flux as a part of energy flux, which is reasonable 
from the relativistic point of view. The energy-momentum tensor can be symmetric 
as we have examined with an MHD wave here. 

What we have shown in the present letter is that the Minkowski momentum can 
be self consistent description of the total momentum of an electromagnetic wave in 
a polarizable medium. This does not necessarily mean the Abraham momentum is 
wrong and inconsistent; it might be possible to give Abraham momentum another 
appropriate meaning to make it consistent. For example, Bamett lUl recently argued 
both Abraham and Minkowski momentum can be consistent when we interpret the 
former as kinetic momentum and latter as canonical momentum. It is out of our 
scope here to examine this argument, however, it should be noted the legitimacy 
of the Minkowski momentum does not automatically exclude the validity of the 
Abraham momentum. 



Appendix 



This appendix is to derive (|T0|) and ((12)) in a shortest way for a physicist not familiar 



with plasma physics. For further information, see any textbook on plasma physics, 
G-gv im Sol- Note that many books derive Alfven waves from the MHD equations, 
which is different from the derivation here; of course the result is the same. 

Suppose a plasma consists of equal number of protons and electrons in a uniform 
magnetic field, which is in the z direction of Cartesian coordinates. The plasma 
response to the electromagnetic field can be expressed by P only and we do not need 
the magnetization current for our calculation. Therefore we can set H = /Jq^B here. 
We assume an MHD wave described above is propagating in this plasma. 

Let us denote a vector in the xy plane by a complex number as A = + iAy. 
Then the wave electric filed in the x direction is denoted as 

E{t) = y (e-'"^ + . (15) 

The equation of motion of a plasma particle in the xy plane is written as 

di? e 

= iQv + -E{t) , (16) 
at m 

where Q. = eB /m is the gyro frequency. We include the sign of the charge in Q, thus Q 
is positive/negative for a proton/electron. The above equation can be directly solved 
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as 



(17) 

2m\D.-co D. + CO 



where Vq is the integration constant. 

Now we assume the wave frequency is much smaller than the gyro frequency 
{co « Q), which is true for most of MHD waves. Then the effect of the first term in 
(fTT)) will be averaged out for MHD time scale; we do not pay attention to this term 
hereafter. 

The rest of the motion is called "drift" in plasma physics. The drift velocity Vd 
can be expanded as 

Vd = ^(icoscot + ^sincot -\ \. (18) 



B \ Q 

The first term of (|T8l) is called E x B drift; protons and electrons drift in the same 
direction with the same speed with this drift. This term is pure imaginary, which 
means the drift is in the y direction. This drift gives the predominant motion of the 
bulk plasma as in ((121), however, it does not cause a current because both protons 
and electrons have the same drift velocity. What contribute to a current is the second 
term of (|T8)) , which is called the polarization drift. Since this term contains Q"^ 
factor, protons and electrons moves in the opposite direction with different speed. 
The electron gyro frequency is much larger than that of protons, therefore, protons 
predominantly carry currents. The drift direction is the same as the electric field since 
it is pure real, and the drift speed is proportional to time derivative of the field because 
of the factor oj and cos cot — > sin cot. Multiplying the proton's second term of (|T8|) with 
the number density and charge, and replacing the factor co and cos cot — > sin cu^ by the 
time derivative, we obtain (|T8)) . 

From the Maxwell's equation we have 

V X V X E = c-^d^E/d^t + ^oj • (19) 

When we assume the wave propagation is in the z direction {d/dx = d/dy = 0) and 
use (|T0|) , we obtain the propagation equation of an MHD wave (Alfven wave) as 

/I l\d^ 

The Alfven speed Va is often much smaller than the speed of light in space and 
laboratory plasmas. We obtain an wave propagating with the Alfven speed when 
we ignore the 1 /c^ term in the above expression. 
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